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Abstract 

Dispersion relations provide a powerful tool to describe the 
electromagnetic form factors of the nucleon both in the 
spacelike and timelike regions with constraints from uni- 
tarity and perturbative QCD. We give a brief introduction 
into dispersion theory for nucleon form factors and present 
results from a recent form factor analysis. Particular em- 
phasis is given to the form factors in the timelike region. 
Furthermore, some recent results for the spacelike form 
factors at low momentum transfer from a ChPT calculation 
by Kubis and MeiBner are discussed. 

1 INTRODUCTION 

A detailed understanding of the electromagnetic form fac- 
tors of the nucleon is important for revealing aspects of 
both perturbative and nonperturbative nucleon structure. 
The form factors also contain important information on nu- 
cleon radii and vector meson coupling constants. More- 
over, the form factors are an important ingredient in a wide 
range of experiments such as Lamb shift measurements [[j]] 
and measurements of the strangeness content of the nu- 
cleon [^]. With the advent of the new continuous beam 
electron accelerators like CEBAF (Jefferson Lab.), MIT- 
Bates, and MAMI (Mainz), a wealth of precise data for 
spacelike momentum transfers has become available. Due 
to the difficulty of the experiments, the timelike form fac- 
tors are less well known. While there is a fair amount of 
information on the proton timelike form factors, only one 
measurement of the neutron form factor from the pioneer- 
ing FENICE experiment exists. 

The basic theoretical framework of dispersion theory for 
nucleon form factors was established long ago [ffl^]. In re- 
cent years, these methods have been refined considerably. 
In particular, constraints from perturbative QCD and inde- 
pendent low-energy experiments (e.g. on the neutron ra- 
dius) have been incorporated [|J [7|] . 

In the first part of this paper, we review the status of dis- 
persion theory for nucleon form factors and show results 
from a recent update In the second part, we dis- 

cuss the effective field theory description of nucleon form 
factors in the framework of ChPT, based on a recent calcu- 
lation by Kubis and MeiBnerfl. 

2 FORMALISM 

Using Lorentz and gauge invariance, the nucleon matrix 
element of the electromagnetic (em) current operator j° m 



can be parametrized in terms of two form factors, 



(p'\jT\p)=u(p') 



F 2 (t) 



u(p), (1) 



where M is the nucleon mass and t = (p 1 — p) 2 = q 2 < 
the four-momentum transfer. F\ and F% are the Dirac and 
Pauli form factors, respectively. They are normalized at 

t = as 

Ff (0) = 1 , F[\Q) = , if (0) = Kp , F"(Q) = n n , 

(2) 

with k p — 1.79 and k„ = —1.91 the anomalous magnetic 
moments of protons and neutrons in nuclear magnetons, re- 
spectively. It is convenient to work in the isospin basis and 
to decompose the form factors into isoscalar and isovector 
parts, 



F? 



1/3 % 



K) , (3) 



where i — 1,2. The experimental data are usually given 
for the Sachs form factors 



G E (t) 
G M (t) 



Fi(t) 
Fi(t) 



rF 2 (t), 
F 2 (t), 



(4) 



where r = -i/(4M 2 ). In the Breit frame, G E and G M 
may be interpreted as the Fourier transforms of the charge 
and magnetization distributions, respectively. 

3 NUCLEON FORM FACTORS IN 
DISPERSION THEORY 

3. 1 Spectral Decomposition 

Based on unitarity and analyticity, dispersion relations re- 
late the real and imaginary parts of the em nucleon form 
factors. Let F(t) be a generic symbol for any one of the 
four nucleon form factors. We write down an unsubtracted 
dispersion relation of the form 

7T J t t'-t-ie. 

where to is the threshold of the lowest cut of F. Since the 
normalization of F is known, a once subtracted dispersion 
relation could be used as well. Eq. ^ relates the em struc- 
ture of the nucleon to its absorptive behavior. The imagi- 
nary part entering Eq. (p|) can be obtained from a spectral 



decomposition \m. For this purpose it is convenient to con- 
sider the em current matrix element in the timelike region, 



Ju 



(N(p)N(p)\f™(0)\0) 
u(p) 



(6) 



Fo(t) 



v{p) 



where p , p are the momenta of the nucleon-antinucleon 
pair created by the current j™ 1 . The four-momentum trans- 
fer in the time like region is t = (p + p) 2 > 0. Using the 
LSZ formalism, the imaginary part of the form factors is 
obtained by inserting a complete set of intermediate states 
as[| 

ImJ M = |(27r) 3 / 2 AA^(p|J,v(0)|A) (7) 

A 

x(X\f^(0)\0)v(p)S i (p + p-p^, 

where TV is a nucleon spinor normalization factor, Z is 
the nucleon wave function renormalization, and Jn(x) — 
JT(x)7o with Jjv(x) a nucleon source. The states |A) are 
asymptotic states of momentum p\ which are stable with 
respect to the strong interaction. For the matrix elements in 
Eq. (^|) to be nonvanishing, they must carry the same quan- 
tum numbers as as the current j^ m : I G (J PC ) = 0~(1 __ ) 
for the isoscalar component and I G (J PC ) = 1 + (1~~) for 
the isovector component of j° m . Furthermore, they have no 
net baryon number. For the isoscalar part the lowest mass 
states are: 3ir, 5ir, . . .; for the isovector part they are: 2tt, 
Air, .... Because of G-parity, states with an odd number of 
pions only contribute to the isoscalar part, while states with 
an even number contribute to the isovector part. Associated 
with each intermediate state is a cut starting at the corre- 
sponding threshold in t and running to infinity. As a con- 
sequence, the spectral function Im F(t) is different from 
zero along the cut from to to oo with to = 4 (9) m 2 for 
the isovector (isoscalar) case. Using Eqs. (^,0), the spec- 
tral functions for the form factors can in principle be con- 
structed from experimental data. In practice, this proves a 
formidable task and is also unstable. However, the spec- 
tral function can be constrained using, for example, vector 
meson dominance. 

3.2 Vector Meson Dominance 

Within the vector meson dominance (VMD) approach, the 
spectral functions are approximated by a few vector meson 
poles, namely the p, ... in the isovector and the u>, (j), . . . in 
the isoscalar channel, respectively. In that case, the form 
factors take the form 



Vi v ' 



(8) 



Clearly, such pole terms contribute to the spectral functions 
as (5-functions, 



These terms arise naturally as approximations to vector me- 
son resonances in the continuum of intermediate states like 
nir (n > 2), NN, KK and so on. If the continuum con- 
tributions are strongly peaked near the vector meson reso- 
nances, Eq. (Q) is a good approximation. 

For the contribution of the two-pion continuum to the 
isovector spectral functions, however, the replacement by 
a sharp p-resonance is not justified and strongly underesti- 
mates the isovector radius. The unitarity relation of Frazer 
and Fulco [ |T(j| ] determines the isovector spectral functions 
from t = 4 m 2 to t ~ 50 m 2 ~ 1 GeV 2 in terms of the 
pion form factor (t) and the P-wave ttttNN partial wave 
amplitudes. The isovector spectral function is strongly en- 
hanced close to the two-pion threshold (see, e.g., Fig. 2 in 
Ref. [^]). The reason for this behavior is well known. The 
corresponding partial waves have a branch point singular- 
ity on the second sheet (from the projection of the nucleon 
pole terms) located at t c — 4m 2 — m*/M 2 « 3.98 m 2 , 
very close to the physical threshold at t = 4m 2 . The 
isovector form factors inherit this singularity and the close- 
ness to the physical threshold leads to the pronounced en- 
hancement. Note that in the VMD approach this spectral 
function is given by a <5-function peak at m 2 ~ 30 m 2 and 
thus the isovector radii are strongly underestimated if one 
neglects the unitarity correction qTT[], as can be seen from 
the formula 



„2\u 



4m; 



dt 
J 2 



(10) 



Im F?' (t) = n a*' S(t~m 2 Vi ) 



Vi 



(9) 



In the isoscalar channel, it is believed that the pertinent 
spectral functions rise smoothly from the three-pion thresh- 
old to the w-peak, i.e. that there is no pronounced effect 
from the three-pion cut on the left wing of the w-resonance 
(which also has a much smaller width than the p). In chiral 
perturbation theory an investigation of the spectral func- 
tions to two loops has indeed shown no such enhancement 
Jl2|], in contrast to the one loop calculation of the isovec- 
tor nucleon form factors where the unitarity correction on 
the left wing of the p has been seen. A form factor anal- 
ysis in the described framework with three vector meson 
poles in both the isoscalar and isovector channels as well 
as the two-pion continuum was carried out by Hohler et al. 

In the following, we describe an extended update that 
also includes timelike data and additional contstraints on 
the spectral function. 

3.3 Constraints 

The first set of constraints concerns the low-i behavior of 
the form factors. First, we enforce the correct normaliza- 
tion of the form factors, which is given in Eq. (^|). Second, 
we constrain the neutron radius from a low-energy neutron- 
atom scattering experiment [pj^.f] 

Perturbative QCD (pQCD) constrains the behavior of 
the nucleon em form factors for large momentum transfer. 

'There has been some controversy about this value recently and for 
future analyses the error of this constraint should be enlarged. 



Brodsky and Lepage [ |l4| | find for t 



Fi(t) - (-t) 



-(*+i) 



111 'a 



1,2, (11) 



where Qq ~ Aqcd- The anomalous dimension 7 depends 
weakly on the number of flavors, 7 = 2.148, 2.160, 2.173 
for Nf = 3, 4, 5, in order. The power behavior of the form 
factors at large t can be easily understood from perturba- 
tive gluon exchange. In order to distribute the momentum 
transfer from the virtual photon to all three quarks in the 
nucleon, at least two massless gluons have to be exchanged. 
Since each of the gluons has a propagator ~ 1/t, the form 
factor has to fall off as 1/t 2 . In the case of F2, there is 
additional suppression by 1/t since a quark spin has to be 
flipped. 

In the form factor analysis, we used spectral functions 
that lead exactly to the large-t behavior as given in Eq. (|ll|). 
The spectral functions of the pertinent form factors are sep- 
arated into a hadronic (meson pole) and a quark (pQCD) 
component as follows, 



F?(t) = F?(t)L(t) 



(12) 



,V. r-l 



L-\m>) 
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A 2 -t 
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Vv 



m v ~ * 



L(t). 



where F[(t) — F[(t) L(i) -1 parametrizes the two-pion 
contribution (including the p) in terms of the pion form 
factor and the P-wave ttttNN partial wave amplitudes in a 
parameter-free manner (for details see [[|). Furthermore, 
the parameter A separates the hadronic from the quark 
contributions. The fits performed in ^ 0] are rather in- 
sensitive to the presence of the logarithmic factor in the 
spectral functions. The factor L(t) in Eq. ( |l2| ) contributes 
to the spectral functions for t > A 2 and in some sense 
parametrizes the intermediate states in the pQCD regime. 
The particular logarithmic form has been chosen for conve- 
nience. More important is the power behavior of the form 
factors, which leads to superconvergence relations of the 
form 



POO 

/ lmFi t2 (t)dt = 0,. 

J to 



(13) 



The asymptotic behavior of Eq. ( Jl 1| ) is obtained by choos- 
ing the residues of the vector meson pole terms such that 
the leading terms in the 1 /i-expansion cancel. 

The number of isoscalar and isovector poles in Eq. ( p"2j ) 
is determined by the stability criterion discussed in detail 
in [^j . In short, we take the minimum number of poles nec- 
essary to fit the data. Specifically, we have three isoscalar 
and three isovector poles. Our fit includes all measured 
form factor data in the space- and timelike regions. We 
stress that we are keeping the number of meson poles fixed 
in order not to wash out the predictive power. Due to the 



various constraints (unitarity, normalizations, superconver- 
gence relations), we end up with only three free parameters 
since two (three) of the isovector (isoscalar) masses can be 
identified with the masses of physical particles. We also 
have performed fits with more poles; these will not be be 
mentioned and we refer the reader to Refs. [^J 0]. 

3.4 Timelike Data 

We should also comment on the extraction of the time- 
like form factor data from experiment. At the nucleon- 
antinucleon threshold, one has by definition 



G M (4Af 2 ) = G£(4M 2 ) 



(14) 



while at large momentum transfer one expects the magnetic 
form factor to dominate. The form factors are complex in 
the timelike region, since several physical thresholds are 
open. Separating \Gm \ and \Ge\ unambiguously from the 
data requires a measurement of the angular distribution, 
which is difficult. In most experiments, it has been assumed 
that either \Gm\ — \Ge \ or \Ge\ — instead. Most recent 
data have been presented for the magnetic form factors. 

3.5 Fit Results 

We discuss two of the various fits performed in Refs. [^, ^] 
in detail: fit 1 is a fit to all spacelike data, while fit 2 also 
includes the data in the timelike region. We stress that in all 
fits, the four form factors have been fitted simultaneously. 
This procedure ensures that the results for different form 
factors are compatible and allows us to test the consistency 
of different data sets. For fit 1 , the vector meson poles in the 
isoscalar channel are the a; (782), 0(1020), and w(1600). In 
the isovector channel, we have the p(1450) and p(1690). 
The third isovector pole is tightly fixed by the constraints; 
its mass turns out to be close to the p(1690). In fit 2, 
all masses are the same, except for the isovector pole at 
y/t = 1690 MeV which is shifted to Vt = 1850 MeV. This 
choice is motivated by the FENICE data which favor an 
isovector resonance slightly below the nucleon-antinucleon 
threshold.0 A detailed listing of all fit parameters and the 
data points included can be found in Refs. [^, M|. 

In Fig. [l], we compare both fits to the world data in the 
spacelike region. The dashed lines indicate fit 1, the solid 
lines indicate fit 2. Both fits give a good description of 
the world data in the spacelike region, differing only for 
momentum transfers where the form factors are not con- 
strained by the data. From fit 1, we extract the following 
nucleon radii [[| 



r p E = 0.847 fm, r 



0.853 fm, 



M 



0.889 fm . 



(15) 



2 Note, however, that despite the fact that the vector meson poles can 
be identified with physical vector mesons, it is not clear whether this in- 
terpretation holds for the higher mass poles. If one tries to fit the masses 
of these poles, they are not well constrained by the data. In any case, these 
poles are a convenient parametrization of the spectral strength in the high 
mass region. 
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Figure 1: Normalized nucleon form factors for spacelike momentum transfer (Q 2 — — t > 0). G 7 ^, G P E , and G V M are 
normalized to the dipole fit, while G E is normalized to the Platchkov fit [ |l5| ] adjusted to give the correct radius. Dashed 
lines: fit 1, spacelike data only. Solid lines: fit 2, including timelike data. 



From the residues of the two lowest isoscalar poles, we can 
extract the ujNN and tfiNN coupling constants [|| R], 



n 2 

y^NN 

47T 

Ait 



34.6 ±0. 



k u = -0.16 ±0.01, (16) 



6.7 ±0.3, k = -0.22 ±0.04, 



where ny is the tensor-to-vector coupling ratio. Both the 
results for the radii and the vector meson couplings are sim- 
ilar to the earlier analysis of Ref . |^] . 

We now turn to the timelike region. In Fig. || we show 
the results of fit 2 compared to the world data. Note that 
all neutron data are from the FENICE experiment [[|] . The 
solid lines show |Gm|> while dashed lines show \Ge\- All 
proton points are for \Gm\ as discussed above. The pro- 
ton magnetic form factor is well described by the fit. The 



FENICE data for the neutron magnetic form factor have 
been analyzed under both the assumption \Ge\ = \Gm\ 
(squares) and \Ge\ = (diamonds). The latter hypoth- 
esis is favored by the measured angular distributions 
Neither data set can be described by the fit. Even when 
the experimental uncertainties of the FENICE data are re- 
duced by a factor of 10 or more, the fit cannot be forced 
through the data. It is possible that the spectral function 
simply has not enough freedom to account for the FENICE 
data and additional poles have to be introduced. It is also 
surprising that the neutron form factor is larger in magni- 
tude than the proton one, as pQCD predicts asymptotically 
equal magnitudes. In any case, there is interesting physics 
in the timelike neutron form factors and more data is called 
for. Finally, it is interesting to note that fit 2 predicts a zero 
in the neutron electric form factor close to threshold. 
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Figure 2: Fit 2 to both spacelike and timelike data for 
timelike momentum transfer. Vertical dotted line indicates 
nucleon-antinucleon threshold. Solid lines show |Gj\/|, 
dashed lines show \Ge\- All neutron data are from the 
FENICE experiment [||. 



4 NUCLEON FORM FACTORS IN CHPT 

Chiral perturbation theory (ChPT) is a low-energy effective 
field theory of the Standard Model. ChPT focuses on the 
universal, low-energy aspects of the physical system. All 
sensitivity to short-distance physiscs is captured in a small 
number of low-energy constants. ChPT allows for a sys- 
tematic and model-independent calculation of low-energy 
observables by means of an expansion in small momenta 
and quark masses. Particularly important is the sponta- 
neously broken chiral symmetry of QCD in the limit of 
vanishing quark masses, which guarantees that this expan- 
sion is well behaved and severely constrains the form of the 
effective Lagrangian. 

The first systematic investigation of nucleon form factors 
in the relativistic version of ChPT was given in Ref. [|T6||. 
This approach, however, has problems with the expansion 
in small momenta since the momenta of the nucleons al- 
ways contain the large nucleon mass. If one uses a nonrel- 
ativistic formulation for heavy nucleons, this problem can 
be circumvented. This Heavy Baryon ChPT is generally 
very successful for one-nucleon observables. The nucleon 
form factors were studied within HBChPT in Refs. [Fl7|. It 
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Figure 3: G E in relativistic ChPT (solid lines) to third 
(upper curve) and fourth order (lower curve). For com- 
parison: heavy baryon approach (lower/upper dashed line: 
third/fourth order), dispersion analysis (dot-dashed curve), 
and some recent experimental points [ 19 1. 



was found that the chiral description already breaks down 
at relatively small momentum transfers of —t = Q 2 ^ 0.2 
GeV 2 . Furthermore, the nonrelativistic expansion does not 
reproduce the correct singularity structure of the isovector 
spectral functions. This problem can be overcome in the 
recently proposed "infrared regularization" of Ref. [|l8|l, 
which is manifestly Lorentz-invariant. Since it is relativis- 
tic, this approach keeps the correct singularity structure in 
the low-energy domain and is expected to improve con- 
vergence of the chiral expansion as well. Recently, this 
method was applied to the electromagnetic form factors of 
the nucleon by Kubis and MeiBner (KM) [§]. In the re- 
mainder of this section, we briefly review their results. 

Fig. |3] shows the electric form factor of the neutron ob- 
tained by KM in relativistic ChPT using infrared regular- 
ization. For comparison, the heavy baryon result and the 
dispersion analysis discussed above (which can be taken as 
a representation of the world data) are also shown. The data 
points are from recent polarization experiments [^9], The 
relativistic ChPT clearly improves over the heavy baryon 
calculation, which already breaks down around Q 2 = 0.1 
GeV 2 . The resummation of the 1/M terms in the relativis- 
tic approach includes recoil corrections and considerably 
improves the convergence of the expansion [j^] . 

In the case of the large dipole like form factors G\, [7 G P E , 
and G P M , however, the relativistic method does not give 
an improvement over the heavy baryon expansion. These 
form factors contain large curvature terms that are not re- 
produced up to fourth order in the chiral expansion. It is 
well known that vector mesons contribute significantly to 
the nucleon form factors and generate the missing curva- 
ture terms. In ChPT these contributions are captured in the 
low-energy constants. The curvature terms, however, only 
appear in higher orders in the chiral expansion. Since the 
vector meson couplings to the nucleon are known from the 
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Figure 4: G\ ;i normalized to the dipole fit in relativistic 
ChPT including vector mesons (solid lines) to third (lower 
curve) and fourth order (upper curve). For comparison: 
world data (cf. Refs. |Q ^) and dispersion analysis (dot- 
dashed curve). 



the form factors over the whole range of momentum trans- 
fers in both the spacelike and timelike regions, while ChPT 
gives a model independent description at low spacelike t 
that explicitly incorporates chiral symmetry. The work of 
KM shows the promise of extending the chiral description 
to higher momentum transfers by merging dispersive meth- 
ods and ChPT [g . 

The neutron form factors in the timelike region show 
some interesting features [ph. In particular, the magni- 
tude of the neutron form factors close to the nucleon- 
antinucleon threshold is larger than expected and at present 
cannot be described by dispersion theory. More precise 
data for the em form factors of the nucleon in the timelike 
region is called for. 
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dispersion analysis described above [|| f7p, it is both reason- 
able and economical to include dynamical vector mesons 
in the theory as no unknown low-energy constants enter. 
Following this spirit, KM include the p(770), w(782), and 
0(1020) as dynamical fields using the technology of anti- 
symmtric tensor fields (see Ref. [|[ for details). Using the 
mechanism of resonance saturation, the low-energy con- 
stants in the theory without vector mesons Li are split into a 
vector meson contribution plus a small remainder Li which 
has to be refit, 



Li 



Li 



(17) 



Keeping the full t-dependence in the vector meson contri- 
bution amounts to resumming a selected class of higher 
order diagrams in the chiral expansion. Fig. |] shows the 
results for G\j obtained by KM using this method. Their 
curves compare well with the world data and the dispersion 
result. KM find that the validity of the chiral expansion can 
be extended up to Q 2 0.4 GeV 2 using this partial resum- 
mation. Similar results hold for the other two dipole-like 
form factors G P M and G P E . Furthermore, the good agree- 
ment for G% is not destroyed by the inclusion of vector 
mesons. 

To summarize, Kubis and MeiBner have shown that one 
can obtain an accurate description of all four nucleon form 
factors up to Q 2 w 0.4 GeV 2 within the framework of rela- 
tivistic ChPT and dynamical vector mesons [H], An exten- 
sion of this work to SU(3) can be found in Ref. [^c||. 

5 CONCLUSIONS 

Dispersion theory and ChPT are two successful approaches 
to the em form factors of the nucleon that mutually comple- 
ment each other. Dispersion theory consistently describes 
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